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The increasing accuracy of experimental data in flavour physics requires a corresponding im- 
provement on the theoretical side, in particular concerning the non-perturbative dynamics of 
QCD. This has prompted the lattice community to aim at an unprecedented accuracy in form 
factors and matrix elements. However, in the light sector, the meson masses remain too heavy 
for an interpolation, which makes it necessary to rely on Chiral Perturbation Theory to perform 
extrapolations in the light quark masses. This makes it all the more necessary to assess precisely 
the range of validity of this theory. More precisely, the presence of strange quark pairs in the 
sea may have a significant impact of the pattern of chiral symmetry breaking : in particular large 
differences can occur between the chiral limits of two and three massless flavours (i.e., whether 
m s is kept at its physical value or sent to zero). We recall some indications of such a scenario 
in QCD, in relation with the peculiar dynamics of the scalar sector. We explain how this could 
affect the convergence of three-flavour chiral series, commonly used to extrapolate the results of 
lattice simulations. Finally, we indicate how lattice simulations with three dynamical flavours 
could unveil such an effect through the quark-mass dependence of light meson masses and decay 
constants. 
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1. Two chiral limits of interest 



Because of the mass hierarchy among light quarks, the strange quark may play a special role 
in the low-energy dynamics of QCD. It is light enough to allow for a combined expansion of 
observables in powers of m u ,md,m s around the Nf = 3 chiral limit (meaning 3 massless flavours): 
m u = ma = m s = 0. But it is sufficiently heavy to induce significant changes in order parameters 
from the Nf = 3 chiral limit to the = 2 chiral limit: m u = rrid = and m s physical. Finally, it is 
too light to suppress efficiently loop effects of massive ss pairs (contrary to heavier quarks) [Eh. 

Two different versions of Chiral Perturbation Theory have been developed around these two 
limits. In the Nf = 2 massless limit, only the pions play a particular role and thus are the only 
available degrees of freedom. The Nf = 3 chiral limit promotes pions, kaons and r\ as the degrees 
of freedom : this second version of #PT is richer, discusses more processes in a larger range of 
energy, but contains more unknown low-energy constants (LECs) and may have a slower conver- 
gence. In each limit, the LECs encode the pattern of chiral symmetry breaking, cannot be computed 
within %?T and must be determined from experiment. Obviously, LECs in both theories are related 
since the Nf = 2 theory can be obtained from Nf = 3 #PT by restricting it to the pion sector and 
integrating out kaons and eta, treated as massive particles. However, the details of the connection 
between the two theories remain under debate. 

Indeed, due to ss sea-pairs, order parameters such as the quark condensate and the pseudoscalar 
decay constant, T.(Nf) = — Mra^ f {uu) and F 2 (Nf) = lim^.F^, can reach significantly different 
values in the two chiral limits (lim#, denoting the chiral limit with Nf massless flavours) An 
illustration is provided by the quark condensate in the two limits: 

1(2) = E(2;m,) =E(2;0) +m / E f ^ + 0{m 2 ) (1.1) 

dm s 

= I(3)+m, lim i f d 4 x (0\uu(x) ss(0)\0) + 0(m 2 s ) (1.2) 

Here, si-pairs are involved through the two-point correlator ((uu)(ss)), which violates the Zweig 
rule in the vacuum (scalar) channel 1 . This loop effect disappears in the large-A^ limit, which is 
known to fail in the ++ channel due to the complicated structure of broad resonances, correspond- 
ing to poles of the scattering matrix located far away from the real axis (see for instance [Q] for a 
discussion of the lightest scalar resonances). 

Arguments based on the spectrum of the Dirac operator [jl|] indicate that this effect should 
suppress order parameters when m s — > 0: E(2) > E(3) and F 2 {2) > F 2 (3). Since the quark con- 
densate^) and the pseudoscalar decay constant(s) are the building blocks of the two versions of 
#PT at leading order, a strong decrease from Nf = 2 to Nf = 3 should have a direct impact on the 
structure of the two theories. We discuss some available data on the pattern of chiral symmetry 
breaking in both chiral limits in turn. 



'Like all the terms in the ntj-expansion of E(2), this correlator is a SUi(2) x SUr(2) order parameter related to the 
spontaneous breakdown of Nf = 2 chiral symmetry, and this for any value of m s . The (scheme-dependent) high-energy 



contributions to eq. (1.2) should be proportional to m = m u = and thus drop from the relation once the Nf = 2 chiral 
limit has been taken (in contrast with the quark condensates arising in OPE at non-vanishing m u ,m^,m s , which exhibit 
such ultraviolet divergences). 
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2. The situation in the Nf = 2 limit 

A few years ago, the Brookhaven E865 collaboration provided new data on decays [Q]. 
Building upon the dispersive analysis of %% scattering [@], we extracted the two-flavour order 
parameters [§]: 

X(2) = h±gp) =0 .81±0.07 Z(2) = ^P= 0.89 ±0.03 (2.1) 

A different analysis, with an additional theoretical input from the scalar radius of the pion, led to a 
larger value for X(2) [g]. 

The situation is somewhat modified by new data from the NA48 collaboration [^], which show 
some discrepancy with the E865 phase shifts in the higher end of the allowed phase space. The role 
of isospin breaking corrections is under discussion currently. The preliminary values of the phase 
shifts [0] tend to increase the value of the / = / = nn scattering length, and to decrease the value 
of the two-flavour quark condensate, possibly pushing X(2) down to 0.7. 

The Gell-Mann-Oakes-Renner relation and its siblings correspond to X(2) =Z(2) = 1 and 
such deviations from unity may seem fairly unimpressive. However one should remember these 
are Nf = 2 chiral expansions in powers of m u and only. Indeed, X(2) and Z(2) indicate the 
convergence of Nf = 2 chiral expansions of F 2 M\ and F 2 respectively : they measure the relative 
size of the leading-order term in these expansions. One usually expects a far quicker convergence 
of Nf = 2 chiral series, with much smaller next-to-leading order corrections (below 10%) 

3. The situation in the Nf = 3 limit 

To include K- and f] -mesons dynamically, one must use three-flavour #PT, where the ex- 
pansion in the three light quark masses starts around the Nf = 3 vacuum m u ,m^,m s = 0. Here, 
strange sea-quark loops may have a dramatic effect on Nf = 3 chiral expansions. The leading-order 
(LO) term, which depends on the 0{p 2 ) low-energy constants F 2 (3) and E(3), would be damped, 
whereas next-to-leading-order (NLO) corrections could be enhanced, in particular those related to 
the violation of the Zweig rule in the scalar sector. For instance, the Gell-Mann-Oakes-Renner 
relation would not be saturated by its LO term and would receive sizeable numerical contributions 
from terms treated as NLO in the chiral counting scheme. 

Indirect estimates |9|] suggest a very significant effect. At NLO the violation of the Zweig 
rule in the scalar sector is encoded in the two LECs L4 and L^. Dispersive estimates of the correlator 
{(uu)(ss)), related to [Bj, indicate that the quark condensate may drop by a half from Nf = 2 
to Nf = 3, i.e. when m s is sent from its physical value down to zero (such a decrease has been 
observed by the MILC collaboration JT0|]). On the other hand 0, the dispersive study of low- 
energy 71K scattering through Roy-Steiner equations and the naive comparison with NLO Nf = 3 
chiral expansions yields U(M p ) = (0.53 ±0.39) • 10~ 3 . 

Such a value may seem rather small, but one should not forget the m s -enhancement of the 
contributions of L4 and L(, in chiral series. Take for instance 

Fl = F(3) 2 + \6(m s + 2m)B AL 4 + \6mB AL 5 + O(m 2 q ) (3.1) 
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where Bo = — ^^m u ,m d ,m s ^o(uu) / F 2 , and we have put together NLO low-energy constants and 
chiral logarithms : AL 5 = L 5 (M p ) + 0.67 • 10~ 3 and AL 4 = L 4 (M p ) +0.51 • 10~ 3 (which is enhanced 
by a factor of m s /m). If we assume that the LO contribution is numerically dominant (i.e., F 2 = 
F(3) 2 to a very good approximation), we can perform the following manipulations: 



F(3) 2 _ F 2 (3) Q 2M 2 +M 2 _ M* 

F 2 -F2(3) + 0(m$)- 1 Fl ^ 4 V 2 



1-8 — K J 7 K M A -%-^^L 5 + 0{m 2 q ) = 1-0.51 -0.04 + 0(m 2 q ) 



(3.2) 

where we have used 1/(1+ x) = 1 — x and eq. ( JO| ) at the second step, and the second and third 
terms of the last equality are obtained using L A {M p ) = 0.5 • 10~ 3 and L 5 (M p ) = 1.4 • 10~ 3 respec- 



tively (these values are used for illustrative purposes). Eq. ( |3.2D is clearly in contradiction with the 
original assumption F 2 ~ F(3) 2 . A similar game can be played with and F 2 M 2 . 

We can draw several conclusions from this simple exercise. A small positive value of L A {M p ) 
or of L(,(M p ) is enough to invalidate the usual assumption of a rapid convergence of Nf = 3 chiral 
series (an issue to be remembered when one tries to extract the values of LECs in lattice simula- 
tions). In addition, we may encounter chiral series A = Alo + Anlo +A8A with a good overall 
convergence, i.e., 8A <C 1 but the numerical balance between LO and NLO depends on the rele- 
vance of strange sea-quark loops. The numerical competition between formal LO and NLO makes 
approximations such as 1 /( 1 +x) ~ 1 — x or (m s + m)2?o — M\ rather hazardous, leading to potential 
contradictions and/or slow convergence for some chiral expansions. 

Actually, such difficulties in the convergence are encountered when NNLO computations are 



fitted to experimental data. A good example can be found among the reference fits in ref . [11]: 
for instance, the so-called Fit 10 yields : (M 2 ) f/ , = (M 2 .)^ [0.736 + 0.006 -0.258], corresponding 
to the relative size of LO, NLO and NNLO respectively. If confirmed, such difficulties in the 
convergence of three-flavour chiral expansions should be considered as a very serious problem and 
a source of sizable systematics in lattice results relying on chiral extrapolations on a large range of 
quark masses. 

4. Constraints from %% and %K scatterings 

We would like to cope with the potentially "large" values of L4 and Lg, and the resulting 
numerical competition between formal LO and NLO contributions in chiral series. To do so, we 



have introduced a framework, called Resummed Chiral Perturbation Theory (Re^PT) []12|], where 



we define the appropriate observables to consider and the treatment of their chiral expansion [[12J, 



13, 14] 



1. Consider a subset of observables that are assumed to have a good overall convergence - we 
call them "good observables". They must form a linear space, which we choose to be that 
of connected QCD correlators (of vector/axial currents and their divergences) as functions of 
external momenta, away from any kinematic singularities. This rule selects in particular Fp 
and F 2 M 2 (P = k,K). 

2. Take each observable and write its NLO chiral expansion in terms of the chiral couplings : 
Fq, Bo, Lj. . . 
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3. In theses formulae, reexpress the bare couplings in terms of physical quantities (masses, 
decay constants. . . ) if justified by physics considerations (e.g., position of unitarity cuts) 
The result is called "bare expansion". 

4. In these bare expansions, reexpress 0(p 2 ) and 0(p 4 ) LECs in terms of: 

m s , , 2m£(3) , , F 2 (3) 

and NNLO remainders, using the bare expansions for the masses and decay constants of the 
pseudoscalar mesons n,K,ri. 

This simple recipe provides a resummation of the potentially large effect of the Zweig-rule 



violating couplings L4 and [|13j, |12j, hence the name of "Resummed Chiral Perturbation Theory" 
(Re^PT) given to this framework. Since we want to cope with the possibility of a numerical 
competition between (formal) LO and NLO terms in chiral series, some usual 0{p A ) results are not 
valid any longer : for instance r = m s /m is not fixed by M\jM\ and can vary from 8 to 40, L5 is 
not fixed by the ratio of Fk/F k . 

In this framework, one can derive the NLO amplitudes corresponding to nn and nK scatter- 
ings, which provide information on Nf = 2 and Nf = 3 patterns of chiral symmetry breaking. The 
smallest uncertainties on the chiral expansion are expected to occur in the unphysical region, far 
away from the non-analyticities due to unitarity. One can exploit dispersive relations, such as Roy 
equations [||] and Roy-Steiner equations [Qj, to reconstruct the amplitudes in this unphysical point 
from the phase shifts from threshold up to energies around 1 GeV. Matching the dispersive and 
chiral representations of the amplitude in a frequentist framework provides constraints (in terms 
of confidence Levels) on the main parameters of interest for three-flavour #PT Jl2|]. As an illus- 
tration, fig. |l| shows the situation for r = m s /m and X(3) = 2m£(3) / (F^M 2 .) . A more detailed 
analysis provides the following constraints from the combination of Tin and nK scatterings: 

r>14.8, X(3)<0.83, 7(3) < 1.1, 0.18<Z(3)<1. [68%CL] (4.2) 
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Fpi"2Mpi A 2 z3 = 0.8 




q 



Figure 2: F% /F% (left) and F%M%/(F%M%) (right) as functions of q = mjm s . Solid, long-dashed and dashed 
curves correspond respectively to X(3) = 0.8,0.4,0.2. Thick (thin) lines are drawn for r = 30 (20). We set 
Z(3) = 0.8 and NNLO remainders are neglected. 



Therefore, one can extract only mild constraints on the parameters of interest from experimen- 
tal data : it is far from clear whether the usual description of chiral symmetry breaking, triggered 
by a large quark condensation and independent from strange sea-quark effects, holds or not. 



5. Three-flavour #PT and lattice simulations 

Lattice simulations may suffer from unexpected systematics in chiral extrapolations due to a 
strong m s -dependence of chiral order parameters, and the resulting numerical competition between 
LO and NLO in three-flavour chiral series. But they may also shed some light on this problem, 
since they allow one to vary the quark masses, and thus enhance or suppress Zweig-rule violating 
contributions accordingly JT4]]. We consider a lattice simulation with (2+1) flavours : two flavours 
are set to a common mass m, whereas the third one is kept at the same mass as the physical strange 
quark m s . Each quantity X observed in the physical situation (m,m,m s ) has a lattice counterpart X 
for (m,m,m s ). One can study the variation of Fp and FpMp according to 

m m s 2m£(3) F 2 (3) 

Figure || illustrates how the dependence of hadronic observables on quark masses is related to 
Zweig-rule violating sea-quark effects: if the latter are large [X(3) = 0.2 <X(2), dashed line], the 
curves bend more than in the case of negligible sea-quark loops [X(3) = 0.8 ~ X(2), solid line]. 

Actually, lattice simulations are performed in a finite spatial box, whereas time is sent to 
infinity to single out the ground state. For sufficiently large boxes, the low-energy effective theory 
is identical to #PT, with the same values of the chiral couplings as in an infinite volume. The 
only difference arises in the propagators of the Goldstone modes, so that the chiral series for Fp 
and FpMp change only through tadpole logarithms. One can extend Re^PT to cope with such a 
problem by noting that only a part of finite-volume corrections computed at NLO in #PT has to 



involve the physical masses of the Goldstone bosons [14]. The relative finite- volume corrections 
for these observables are significant for L = 1.5 fm, but much smaller for L = 2.5 fm. In addition, 
the corrections are smaller for F%M% than for the decay constant. By inspection, we see that FpMp 



6 



Surprises and pitfalls in three-flavour chiral extrapolations 



Sebastien Descotes-Genon 



(P = K, K) in large volumes is a quantity for which we manage a good control of finite-volume 
effects, whatever the impact of sea ss pairs on the pattern of Nf = 3 chiral symmetry breaking 
More specifically, the two ratios: 

IFffil 2 tlMl 

R *-qF&i> RK -(g + l)FiMr ^ 

are only mildly affected by finite-volume effects, and their ^-dependence should provide valuable 
information on the importance of strange-quark loops. NNLO remainders blur slightly the picture, 
but they do not prevent the assessment of the effect. 

6. Conclusion 

The presence of massive 55-pairs in the QCD vacuum may induce significant differences in 
the pattern of chiral symmetry breaking between the Nf = 2 and Nf = 3 chiral limits. This effect, 
related to the violation of the Zweig rule in the scalar sector, may destabilise three-flavour chiral 
expansions numerically and may yield large systematics for chiral extrapolations of lattice data 
down to very light quark masses. We have compared the situation in the two- and three-flavour 
chiral limits : the combination of %% and %K data does not favour the usual picture of a large quark 
condensation independent of the number of massless flavours. We indicated how chiral extrapo- 
lations of lattice results with three dynamical quarks could be affected in the case of significant 
differences between the Nf = 2 and Nf = 3 chiral limits. In addition, we suggested a way to test 
this scenario on the lattice, through two dimensionless ratios which prove sensitive to this effect, 
with only a mild impact of finite-volume corrections. 
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